GRB 010222: A Burst Within a Starburst. 
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ABSTRACT 

We present millimeter and submillimeter wavelength observations and 
near-infrared i^-band imaging toward the bright gamma-ray burst GRB 010222. 
Over seven epochs the flux density of the source was constant with an average 
flux density 3.74±0.53 mJy at 350 GHz and 1.05±0.22 mJy at 250 GHz, giving 
a spectral index a = 3.78 ±0.25 (where Foe z/"). We rule out the possibility that 
this emission originated from the burst or its afterglow and we conclude that 
it is due to a dusty, high redshift starburst galaxy (SMM J14522+4301). We 
argue that the host galaxy of GRB 010222 is the most plausible counterpart of 
SMM J14522+4301, based in part on the centimeter detection of the host at the 
expected level. The optical/NIR properties of the host galaxy of GRB 010222 
suggest that it is a blue, sub-L,,, similar to other GRB host galaxies. This 
contrasts with the enormous far-infrared luminosity of this galaxy based on our 
submillimeter detection (LboI ~ 4 x lO^^L©). We suggest that this GRB host 
galaxy has a very high star formation rate, SFR«i600 Mq yr~^, most of which is 
unseen at optical wavelengths. 

Subject headings: gamma rays: bursts - radio continuum: general 
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1. Introduction 

GRB 010222 was detected on 2001 February 22.31 UT by the gamma-ray burst monitor 
on tiie BeppoSAX satellite and localized to a error circle of 2.5' by the Wide Field Camera 
(Piro 2001a,b). A bright optical afterglow candidate was quickly identified by Henden 
(2001) and subsequently confirmed by follow-up optical and radio observations ( [McDowell 



|2001| , Perger fc Frail 2001| ). Absorption lines measured by Jha et al. (2001) yielded a lower 
redshift limit to the GRB of z = 1.477, which they argued to be the redshift of the GRB 
host galaxy based on the unusually strong equivalent widths of the metallic lines. Since 
the burst and its afterglow were bright we initiated GRB target-of-opport unity programs 
at several different telescope facilities ( Frail et al. 2000| ). 



In this paper we report on observations of GRB 010222, which contrary to expectations, 
was detected as a constant millimeter and submillimeter source while the flux density at 
radio and optical wavelengths, showed the usual declining behavior expected from the 
afterglow. After considering several different origins for the emission, we conclude that we 
have detected the host galaxy of the burst, which radiates primarily in the infrared (IR) 
and submillimeter, due to the dust and gas heated by intense star formation. 



2. Radio Observations 

James Clark Maxwell Telescope (JCMTf\): A frequency-ordered log of all observations and 
flux density measurements are given in Table |^. Observations were made in the 350 GHz 
and 660 GHz bands using the Sub-millimeter Common-User Bolometer Array (SCUBA; 
Holland et al. 19991) . In this mode, the secondary mirror is chopped between the on-source 



position and a position 60" away in azimuth, in order to remove sky emission. The angular 
resolution of the telescope is ~14" at 350 GHz and ~8" at 660 GHz. The pointing was 
checked approximately once per hour on a nearby blazar (JVAS 1419+5423) which is 
located less than 13° from the GRB. The pointing was generally found to vary by less than 
2-3". On all days where observations were taken, the focus was checked at least every 2-4 
hours, usually on the (local) pointing source. 



The data were reduced using the SCUBA User Reduction Facility ( [Jenness fc Lightfoot 



19981) . The raw signals were fiat-fielded to account for the small differences in bolometer 



®The JCMT is operated by The Joint Astronomy Centre on behalf of the Particle Physics and Astronomy 
Research Comicil of the UK, the Netherlands Organization for Scientific Research, and the National Research 
Council of Canada. 



- 3 - 



response, extinction corrected, and de-spiked to remove anomalous signals above the 
3-sigma level. Short time-scale sky variations were also removed using pixels around the 
edge of the array containing no source emission ( [Jenness, Lighttoot fc Holland 19"^ ). A 



flux calibration factor was then applied to convert to Jy. Flux calibration factors of 197 ± 
13 Jy/V and 384 ± 82 Jy/V were applied to the 350 GHz and 660 GHz data, respectively 
( IGoulson 200"0|) . 

The atmospheric opacity at 220 GHz was measured every 10 minutes using a tipping 
radiometer mounted near and operated by the Caltech Submillimeter Observatory. These 
values were then extrapolated to 350 GHz and 660 GHz and used to correct for atmospheric 
extinction ( [Archibald, Wagg fc Jenness 2 



IRAM 30-meter telescope (IRAM): Observations were made using the Max-Planck 
Millimeter Bolometer (MAMBO; [Kreysa et al. 1998|) array at the IRAM 30-m telescope on 



Pico Veleta, Spain. MAMBO is a 37-element bolometer array sensitive between 190 and 
315 GHz. The larger than half peak sensitivity range is 210 - 290 GHz, with an effective 
center frequency of ~ 250 GHz for steep spectrum sources. The beam for the feed horn of 
each bolometer is matched to the telescope beam of 10.6", and the bolometers are arranged 
in a hexagonal pattern with a beam separation of 22". Observations were made in standard 
on-off mode, with 2 Hz chopping of the secondary by 32". 

The telescope pointing was checked every ~ 20 minutes with JVAS 1506+4239, a point 
source ~2.5 deg from the OT position. The flux of the pointing source was monitored and 
appeared stable at 0.15 Jy. The target pointing accuracy is better than 2". GRB 010222 was 
positioned on the central bolometer of the array, and after each 10 seconds of integration, 
the telescope was nodded so that the previous "off" beam becomes the "on" beam. Each 
scan of sixteen 10 second subscans lasts 4 minutes, of which 80 sec integration falls on 
the source, 80 sec off source, and about 80 sec are lost to move the telescope and start 
integration. Gain calibration was performed using observations of planets, resulting in a 
flux calibration factor of 12,500 counts per Jansky, which is estimated to be accurate to 
15%. Opacity corrections were made from frequent skydips. A gain-elevation correction 
was applied. 

The data were reduced using the MOPSI software package (Zylka 1998). The average 
flux densities found at the target position along with other observational parameters are 
given in Table |l|. 

Owens Valley Radio Observatory Interferometer (OVRO): Observations were made on 
two nights at a central frequency of 98.48 GHz with a 1 GHz wide band. On the second 
night, weather conditions improved sufficiently to add a second frequency at 220.73 GHz 
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with a 2 GHz wide band. Phase calibration was performed using the source 3C 345, while 
observations of Neptune provided the flux density calibration scale with an estimated 
uncertainty of ~ 15%. See Shepherd et al. (1998), for further details on calibration and 
imaging. No source was detected on either night. 



3. Optical and Near IR Observations 

MDM Observatory: The position of GRB 010222 was observed in the K band using the 
TIFKAM infrared imager/spectrograph on the 2.4 m Hiltner telescope at MDM Observatory 
on 2001 April 11. This instrument uses a 512 x 1024 InSb detector, and covers a 2^5 x 5' 
field of view. A total of 135 dithered 30 sec exposures were sky-subtracted and combined 
to yield a total exposure of 67.5 min. The seeing was 0."8 — 0."9, and the airmass ranged 
from 1.02 — 1.14. Nine standard stars from the list of Persson et al. (1998) were observed 
throughout the night. From these observations we find that the photometric zeropoint 
accuracy is 0.1 mag due, probably, to the presence of light cirrus. 

The magnitudes of several objects in the neighborhood of the GRB host galaxy were 
measured. For reference, a magnitude K = 16.507 ± 0.033 was determined for the star "E" 
(Figure]^) which lies 10" south and 5.3" east of GRB 010222. Based on this calibration and 
the photometry of the faintest visible objects near the OT, a Aa upper limit K > 19.6 was 
derived at the position of GRB 010222. 

Keck Observatory: The location of GRB 010222 was observed with the NIRSPEC 
slit-viewing camera on Keck 11^ on 2001 April 27 UT, approximately 64 days after the GRB. 
The observations consisted of a total of 58 minutes integration over 20 individual dithered 
frames with the airmass ranging from 1.3 - 1.7. The conditions were photometric with 
~ 0.45" seeing. The images were sky-subtracted and combined using the IRAF/DIMSUM0 



The MDM-calibrated magnitude of star E (see above) was used as the basis of the 
magnitude zero-point. The host of GRB 010222 is detected at the ~ 3.6a level (see Figure 
For completeness the four other objects in the i^-band image that lie within a 7" 



^The W. M. Keck Observatory is operated by the California Association for Research in Astronomy, a 
scientific partnership among CaHfornia Institute of Technofogy, the University of Cahfornia and the National 
Aeronautics and Space Administration. 
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radius are indicated in Figure 0. Table contains the results of the i^-band and F814W 
photometry and the (F814W— i^) colors (see below). Aperture photometry was performed 
on each of these sources with circular apertures whose radius was chosen to include most of 
the object flux, based on a curve of growth. 

In order to derive colors for these galaxies matching apertures were used (with 
re-centering allowed) to photometer HST+WFPC2 F814W images of the field. These HST 
data were taken over four epochs and are part of a multi-epoch, broadband monitoring of 
the afterglow of GRB 010222. Full details of the HST observations can be found in Galama 
et al., {in prep). 

4. The Origin of the Submillimeter and Millimeter Emission 

Emission was detected towards GRB 010222 at 250 GHz and 350 GHz (see Table |lD. 
In contrast to the steady decay of the centimeter and optical emission, the submillimeter 
flux at 350 GHz remains remarkably constant from 6 hours after the burst, until 18 days 
later (Figure p. The noticeable exception to this behavior is on Day 8, when the measured 
flux density lies 2.4(T below the average of all measurements. The data were examined 
carefully for errors and inconsistencies in the calibration and data taking, but apart from 
some problems with the telescope focus at this time (R. Ivison, priv. comm.), we can find no 
instrumental origin for the deviation. Nevertheless, significant fiuctuations in photometric 
measurements taken with single dish submillimeter telescopes are not uncommon. For 
example, the 350 GHz photometry of the z = 1.44 ultraluminous starburst galaxy ERO 
J164502+4626.4 reported by Cimatti et al. (1998) is twice the correct value (Dey et al. 
1999). A search for variability was carried out within individual days but no statistically 
significant evidence for any variations was found. Adding the 350 GHz data from all days 
together gives an average fiux of 3.13±0.47 mJy. Excluding the low measurement gives an 
average of 3.74±0.53 mJy. 

Independent support that the source is constant comes from the IRAM 250 GHz 
measurements. Although the fiux level is weaker than at 350 GHz (Figure |1|), the average 
fiux density at 250 GHz over the same time interval (found by averaging all observations 
while weighting each scan with the inverse square of its average noise level) was 1.05±0.22 
mJy. A later measurement in 2001 May 5 gives a consistent value of 1.13±0.48 mJy. The 
350 GHz to 250 GHz fiux ratio of these average fiux density values implies a steep spectral 
index a = 3.78 ± 0.25 (where Fcx z/"). 

The constancy and steep spectral slope of the submillimeter and millimeter emission 
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towards GRB 010222 are anomalous compared to the behavior of past GRBs and when 
compared to observations of this event at other wavelengths. In particular, the temporal 
behavior is difficult to reconcile with the expectations of the simplest afterglow models. 
This point is illustrated in Figure |T] where we show the deviation of the 250 GHz and 350 
GHz measurements from the predicted afterglow light curves, as determined by broadband 
fitting of the radio, optical and X-ray data (Galama et al., in prep). 

Similar difficulties are posed by the peculiar spectral slope (see Figure 0), whose value 
between 250 GHz and 350 GHz is a = 3.78 ± 0.25. Free-free or synchrotron self-absorption 
( |Waxman 1997| , [Katz fc Piran 1997D in the millimeter band can lead to steep slopes with 



a = 2 — 2.5. However, this explanation fails because the source was readily detectable at 
centimeter wavelengths (|Berger fc l^rail 2001| ) where a simple extrapolation of this slope 
should suggest otherwise. 

Such steep slopes are reminiscent of thermal emission from dust. At wavelengths 
longer than optical, the emissivity of dust grains is usually modeled as Q{\) oc where 
1 < /5 < 2. For dust which is optically thin in the Rayleigh- Jeans portion, this results is a 
spectrum of F^, oc z/^+^, consistent with the observed steep slope. We begin by considering 
whether the dust could have been heated (and partially evaporated) by the burst itself or 
by the prompt optical/UV flash generated in the reverse shock of the afterglow (fWaxman 



|fc Draine 2000| , [Esin fc Blandford 2000| , [Reichart 200l| , [Venemans fc Blain 200 The 



subsequent reprocessing of this radiation, either from dust scattering or by thermal emission, 
results in a long-lived infrared emission component (aka dust echo). A simple energy 
calculation demonstrates the difficulties with this interpretation. For our adopted cosmology 
(Ho=65 km s"^ Mpc~^, fiM=0.3, and Ao=0.7) the luminosity distance corresponding to the 
lower redshift limit of z = 1.477 ( |Jha et al. 2001| ) is di = 3.56 x 10^^ cm. The isotropic 
luminosity at 350 GHz oc u x F,^ (defined in the rest frame) is Lj^ ^ 5 x 10^^ erg s~^. 
Within 8 days of observation, the total energy emitted directly into the 350 GHz band is 
of order 10^° erg. This dust emission is expected to be roughly isotropic, so there should 
be no geometric correction even if the burst itself was beamed. More importantly, the 
bolometric corrections are expected to be huge. Calculations by Waxman & Draine (2000) 
and Venemans & Blain (2001) show that dust is sublimated by the optical/UV flash at any 
distance closer than about 10 pc, and the surviving dust just outside this distance is heated 
to 2000 K, with a spectral peak at a rest-frame wavelength of 1-10 /im. In this case the 
bolometric correction to L^, at 350 GHz is more than six orders of magnitude. Even in the 
most optimistic case, the optical/UV flash is only expected to be a small fraction of the 
GRB energy ( ^ari fc Piran 1999| , [Waxman fc Draine 2000| ) and therefore it is not enough to 
account for the total energy radiated at 350 GHz during the duration of our observations. If 
one hypothetically sets up a geometrical situation in which the dust is located farther away 
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from the burst, to allow for lower temperatures and therefore smaller bolometric correction 
to the 350 GHz flux, one finds that the emission will last for years, increasing again the 
required energy. We therefore conclude the GRB is not energetic enough to create the 
sub- mm flux by heating its surrounding dust. 

We are led to a model in which the millimeter and submillimeter emission is not 
dominated by the direct or reprocessed light from the afterglow itself, but rather originates 
as a distinct and separate component. The simplest hypothesis, consistent with the 
persistent nature of the source and its steep spectrum, is that we have detected a dusty, 
high redshift galaxy coincident with GRB 010222, designated as SMM J14522+4301. The 
population of submillimeter galaxies, of which SMM J14522+4301 is an example, produce 
the bulk of the submillimeter background and are an important contributor to star 
formation at high redshifts ([Small, Ivison fc Blain 19"W[), such systems appear to be analogs 



of local ultraluminous infrared galaxies (ULIRGs, e.g., Sanders fc Mirabel 1996| ). In the 
next two sections we will examine this hypothesis in more detail and derive the physical 
properties of galaxy. 



5. Submillimeter and Millimeter Galaxy Properties 

If GRB 010222 occurred in SMM J14522+4301, then its physical properties can be 
derived. The FIR and submillimeter emission from SMM J14522-I-4301 originates from dust, 
heated by an interstellar radiation field dominated by massive stars ( [Kennicut 1998| ). In 
this instance we are well-justified in representing the thermal dust emission spectrum by a 
modified blackbody which takes the form in the Ray leigh- Jeans regime of F,^ oc z/" x By{T). 
A fit of the OVRO, IRAM and JCMT flux density measurements (and upper limits) results 
in a bolometric luminosity of LboI = 4 x 10^^ Lq. A dust temperature of T(i=40 K and 
a dust emissivity index /3 = 1.5 were assumed ( |Blain l9^ ). For a continuous starburst, 
the star formation rate SFR~ (5mf(Lboi/1O"'^°L0) M© yr^-*^ ( Pmont et al. 200ID , where ^mf 



depends on the stellar mass function. At starburst ages of 10-100 Myr, 5uf lies between 0.8 
and 2, so adopting 5mf=1-5 we derive SFR~600 M© yr~^ Although the uncertainties are 
at least a factor of two, these values of SFR and LboI place SMM J14522+4301 in the class 
of ULIRGs (ISmail, Ivison fc Blain 1997|). 



With what degree of confidence can we assume that SMM J14522-I-4301 is the host 
galaxy of GRB 010222? From submillimeter source counts of dusty galaxies (p31ain et al. 



|1999|) background source confusion is expected at the flux level of SMM J14522+4301 
in the 14.5" (FWHM) JCMT beam at 350 GHz 2±1% of the time. Given the surfeit 
of NIR/optical candidates in Figure ^ it is not sufficient to secure the identification of 
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GRB 010222 with SMM J14522+4301 on a posteriori grounds alone. Similar claims have 
been made for an ISO sateUite detection of a persistent source towards GRB 970508 at 60 
fim (|Hanlon et al. 20001). In our view, independent corroborative evidence is needed before 



any such associations can be made. 

The problem here is a common one. Because of the large beams used in their detection, 
source confusion poses considerable difficulties for identifying the optical counterparts 
of submillimeter-selected galaxies [e.g., [Small et al. 20001) . In addition to a detection 



coincident with GRB 010222, there are four other objects in the K-hand image that lie 
within the JCMT beam (Figure ^ and Table 0). Deep observations (/ > 26 and K > 21) 
towards submillimeter sources have revealed a growing number of highly obscured sources 
(EROs) with very red colors (/ — A") > 5 ( [Smail et al. 199% pey et al. 199% [brayer et~al. 



120001 ). Object "B", 4" away from GRB 010222 is the reddest galaxy in the error circle with 
(F814W— -A) ^ 4.5, nominally qualifying it as an EROs. The red colors of ERO galaxies 
may be due to an old stellar population or a highly reddened population of young stars. 
Their space density, although not well determined, are comparable to that of SCUBA 
sources at these flux levels ( |Uowie et al. 199^ , Moustakas et al. 1997| , Small et al. 1999|) . 



However, not all submillimeter galaxies are red objects. A case in point is 
SMM J14011+0252 at ^ = 2.56, ( [Ivison et al. 2000|) which consists of two components 
(Jl and J2) neither of which are especially red (/ — K) ~ 3. Still others, such as 
SMMJ02399-0136 at z = 2.80 ([Ivison et al. 1998|) , and SMM J02399-0134 at ^ = 1.06 
( poucail et al. 19"99|) have blue colors. This diversity of optical properties is also reflected in 
samples of local ULIRGs, which have SEDs that when shifted to high z can be as blue as 
optically-selected field galaxies ( [Trentham, Kormendy fc Sanders 1999[ ). 



The detection of centimeter radiation coincident with submillimeter sources has proven 
to be a powerful method for identifying counterparts ( [Small et al. 2000| ). Likewise, the 
most secure identification of an ULIRG with a GRB is that of Berger et al. (2001), where 
a persistent centimeter source was found spatially coincident with GRB 980703. The 
method relies on the well-known radio/FIR correlation for star-forming galaxies ( [Condon 



1992). The thermal dust emission at 350 GHz due to massive stars can be related to the 
synchrotron and thermal Bremsstrahlung emission at 1.4 GHz via a redshift-dependent 
spectral index af^^ ( [Carilli fc Yun 19"99[ ). Using the relation of Carilli & Yun (2000) we 



find a value of af^^ = 0.69 ± 0.16 at z = 1.477, and predict a centimeter fiux density at 
1.4 GHz of 8515^^ /iJy. Galama et al. (2001) have noted a fiattening of the 8.46 and 4.86 
GHz light curves at late times. This deviation from model predictions of the afterglow fiux 
is significant at the 2.5-0" and 2.3-cr levels, respectively. For a centimeter spectral index of 
a = —0.7, typical of such galaxies, we find that within the errors the centimeter emission is 
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consistent with an ULIRG at z ~ 1.5. There are no other radio sources visible in the JCMT 
beam to an rms ~6 /iJy beam~^ with a ~1" beam at 8.46 GHz. Continued monitoring will 
verify whether this is the centimeter detection of the host galaxy of GRB 010222. 

In the absence of suitable optical/NIR color discriminants, submillimeter sources have 



been identified by searching for molecular gas at the candidate redshift {^e.g., [Frayer el 



al. 1998|), or by better constraining the position with a interferometric detection of the 



millimeter continuum {e.g., [Bertoldi et al. 2000| ). If SMM J14522+4301 is an ULIRG at the 



same redshift of GRB 010222, then it is estimated from the observed flux density at 250 
GHz (and assuming Trf=40 K, P = 1.5, dust-to-gas ratio=100) a total gas mass of 3x10^° 
Mq (see Omont et al. 2001 and references therein) - an amount that is readily detectable 
with current instruments. A future search for redshifted C0(2-l) at z = 1.477 would be a 
promising means to confirm the association. 



6. Optical and NIR Galaxy Properties 

The current evidence identifies the host galaxy of GRB 010222 as the most likely 
counterpart of SMM J14522+4301. The argument is supported on probabilistic grounds 
and more importantly, the tentative detection of the host at centimeter wavelengths 
(§^. The optically visible portion of the host galaxy of GRB 010222 is blue with 
(F814W— fC) 2.1. The rest-frame B-band luminosity of this galaxy [z = 1.477) is 
Mb = -18.6 (or Mab{B) = -18.8). We have used the formalism of Lilly et al. (1995), 
and assume C(B,I)=— 0.2 for the color term of the K-correction, appropriate for late-type 
galaxies. This is similar to the host galaxy of GRB 970228 ( [Bloom, Djorgovski fc Kulkarn^ 
2000| ) both of which are subluminous (L ~ O.IL^.) relative to L^, galaxies at their respective 



redshifts. However, note that neither the magnitude nor the luminosity of this galaxy is 
exceptional compared to other GRB host galaxies ( [Mao fc Mo 1998| , [Hogg fc I'ruchter 



[1999[ , pchaefer 2000[) . The Mb luminosity corresponds to a rest-frame blue luminosity 
Lb — X L,y = 2 X 10^ L0, which is small compared to the bolometric luminosity of 
LboI = 4 X 10^^ L0 derived in Despite the prodigious star formation, the distribution of 
dust in the host galaxy must consist of regions of significant internal extinction while still 
allowing the escape of some UV light. 
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7. Discussion and Conclusions 



Multi-wavelength studies of GRBs and their host galaxies have the potential both to 
determine the fraction of star formation that is optically obscured ( |Wijers et al. llJ95| , 



Totani 1999|, |Blain fc Natarajan 2000|) and to elucidate the relation between the optical/UV 
and submillimeter selected galaxies. To the degree that GRBs trace massive star formation 
in the early universe ([Bloom et al. 19"99| , [Galama et al. 2000| , Piro et al. 2000| , pioom, 



Kulkarni fc Djorgovski 2000 , Reichart 2001 ), then the extreme luminosity and the dust 
penetrating ability of 7-rays can be used to select high redshift galaxies in an unbiased 
way without regard to their emission properties. Surveys for optical/UV galaxies are 
magnitude-limited while, in contrast, the sample of ORB host galaxies have magnitudes 

~ 22 — 28 which span the complete range of model predictions ( [Hogg fc Fruchter 1999| ). 
Obtaining reliable galaxy identifications is a ongoing problem for submillimeter sources but 
the early detection of an afterglow ensures that both an accurate position and (in many 
cases) a redshift will be known for the galaxy. 

In this paper we have identified a persistent, steep spectrum submillimeter source and 
concluded that it is a dusty, starburst galaxy, most likely the host galaxy of ORB 010222. 
Although centimeter and submillimeter observations of GRB host galaxies have only just 
begun ( [Smith et al. 199U[ , [Berger, Kulkarni fc I'rail 200 1[ , [Hanlon et al. 2000[ ), they are 
revealing aspects of the sample that were not apparent from more extensive optical/NIR 
observations. The detection of significant submillimeter and centimeter radiation suggests 
that at least some GRB host galaxies are undergoing prodigious star formation (> SOOM© 
yr~^). Such results are important as they are a direct link between GRBs and episodes of 
massive star formation (Berger et al. 2001). 

Future deep submillimeter and centimeter radio observations of GRB host galaxies 
hold some promise for exploring these differences in their optical/UV properties. The 
current sample is somewhat limited. Centimeter observations have not been carried out 
to the requisite level {i.e., at z=l, 10 yuJy— IOMq yr~^ at 1.4 GHz) and to do so for the 
entire sample would require the sensitivity of the Expanded VLA (EVLA), currently under 
development. Submillimeter and millimeter observations have been undertaken toward 
10 well-localized GRBs ( Shepherd et al. 1998[ , [Bremer et al. 1998[ , [Kulkarni et al. 



1999 



Smith et al. 1999| , [Frail et al. 2000| , [Berger et al. 2000[) but for the most part the existing 
upper limits would not be sufficient to rule out galaxies as bright as SMMJ14522-I-4301. 
Ramirez-Ruiz et al (2001) have calculated that 20% of GRB host galaxies are detectable at 
350 GHz with the SCUBA array on JCMT, a number that is close to the current detection 
rate. 
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Table 1. Submillimeter and Millimeter Observations of GRB 010222' 





Freq. 


Start 


UT Start-End 


TOS 


Fluxibrms 




Notes 




Telescope 


(GHz) 


Date 


(H:M) 


(s) 


(nijy) 








JCMT 


660 


Feb. 


24 


16:17- 


17:56 


3600 


< 37.8 (3cr) 


r 


= 0.055 




JCMT 


350 


Feb. 


22 


13:03 - 


18:39 


10800 


4.34 ± 1.09 


r 


= 0.11 




JCMT 


350 


Feb. 


23 


11:08 - 


19:30 


18900 


3.36 ± 0.72 


T 


= 0.09 




JCMT 


350 


Feb. 


24 


16:17- 


19:29 


6300 


4.04 ± 1.14 


T 


= 0.055 




JCMT 


350 


Mar. 


01 


11:54- 


12:50 


2070 


1.14 ± 2.06 


T 


= 0.08 




JCMT 


350 


Mar. 


02 


14:54- 


18:05 


7200 


0.63 ± 0.94 


T 


= 0.05 




JCMT 


350 


Mar. 


12 


16:18 - 


19:07 


5400 


3.98 ± 1.25 


T 


0.05 




IRAM 


250 


Feb. 


22 


22:33 - 


06:43 


8850 


1.32 ± 0.54 


T 


= 0.45 - 


0.51 


IRAM 


250 


Feb. 


25 


00:28 - 


07:55 


9140 


1.09 ± 0.32 


T 


= 0.17- 


0.29 


IRAM 


250 


Feb. 


26 


01:12 - 


08:52 


8160 


1.00 ± 0.33 


T 


= 0.23- 


0.28 


IRAM 


250 


Mar. 


13 


07:55 - 


08:09 


660 


0.91 ± 0.94 


T 


= 0.08 




IRAM 


250 


Mar. 


14 


05:31 - 


06:00 


1330 


1.31 ± 0.61 


T 


= 0.12 




IRAM 


250 


Mar. 


18 


22:33 - 


22:48 


830 


0.78 ± 1.10 


T 


0.17 




IRAM 


250 


May 05 


22:57 - 


03:18 


2300 


1.13 ± 0.48 


T 


==0.14- 


0.16 


OVRO 


221 


Feb. 


23 


07:26 - 


17:40 


18000 


1.15 ± 5.00 


T 


= 0.25 




OVRO 


98.5 


Feb. 


22 


14:31 - 


19:44 


8400 


0.30 ± 2.00 


T 


> 0.25 




OVRO 


98.5 


Feb. 


23 


07:26 - 


17:40 


18000 


-0.30 ± 0.80 


T 


= 0.25 





^From left to right the columns are the telescope used, the sky frequency at which the observations were 
made, the starting date, the UT start and end time, the time on source in seconds, the peak flux density 
at the position of GRB 010222 and the rms noise in milliJansky, and the zenith atmospheric opacity at 220 
GHz (250 GHz at IRAM). See §^ for more details on the data calibration. 
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Table 2. Magnitudes and Colors of K-band selected galaxies within the SCUBA beam 



Galaxy 


Aperture (") 


F814W 


K" 


F814W-/f 


A 


1.42 


21.786±0.010 


18.233±0.034 


3.553±0.035 


B 


1.60 


22.931±0.028 


18.398±0.035 


4.533±0.044 


C 


1.96 


22.124±0.018 


19.049±0.041 


3.074±0.044 


D 


1.42 


24.934±0.151 


21.102±0.113 


3.832±0.188 


Host 


0.45 


25.58±0.13 


22.745±0.274'= 


2.834±0.303 



'^The error estimates do not include a zeropoint uncertainty of 0.1 mag (see 

^The host galaxy of GRB 010222 in the K band is contaminated to some extent by the afterglow. An 
estimate from broad band modeling is that the afterglow contributes approximately 50% of the total flux 
at this epoch. There is no afterglow component in the F814W band, since the host and afterglow were fit 
independently. After correcting for the afterglow contamination we find that the host color is bluer, i.e. 
(F814W-i4:) ~ 2.1. 
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Fig. 1. — Flux density measurements toward GRB 010222 taken with the SCUBA bolometer 
array on the James Clerk Maxwell Telescope (top panel) and with the Max-Planck Millimeter 
Bolometer array (MAMBO) on the IRAM 30-meter telescope (bottom panel). Upper limits 
are shown as open triangles and are plotted as the peak flux density at the location of the 
afterglow plus two times the rms noise. The observing frequency is shown in the upper left 
corner of each panel. The dotted lines are afterglow hght curves derived from a global fit to 
the entire broad-band dataset (see Galama et al., in prep). 
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Fig. 2. — Spectral energy distribution (SED) towards GRB 010222. The dashed line is 
the contribution of the ULIRG SMM J14522+4301 to the total emission. It consists of a 
modified blackbody of the form F^, oc B,^{T){1 — exp[— (z//z/o)°]} (where T=40 K, a = 1.5 
and z/q = 10^^ Hz) plus a power-law synchrotron spectrum of the form F^, oc i/" where 
ar = —0.75 (see text for more details). To this starburst SED we have added the SED for 
the afterglow emission from GRB 010222 (solid lines) at two epochs, 2.4 and 23.6 days after 
the burst. These have been derived from broadband modeling of all available radio, optical 
and X-ray measurements taken between 0.15 days to 75 days after the burst (Galama et al., 
in prep). The millimeter and submillimeter flux density measurements are from Table |I|, 
the centimeter radio and NlR/optical data are from Galama et al., and the IRAS data were 
obtained from the NASA/IPAC InfraRed Science Archive. To better constrain the SED of 
SMM J14522+4301 we have used the average flux densities at 250 GHz (IRAM) and 350 
GHz (JCMT), rather than single epoch measurements. 
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Fig. 3.— K-band Keck II+NIRPSEC image of the field of GRB 010222. The 7-arcsec radius 
of the SCUBA beam is indicated and the four galaxies (A-D) discussed in the text are labeled. 
Also marked is the 3a error ellipse of the optical transient based on a relative astrometric 
solution with the HST+WFPC2 third epoch F814W observation of the afterglow. The host 
of GRB 010222 is detected at the 3.6a level. North is up and east is to the left. 



